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In  this  paper  we  propose  a  fully  conservative  torm  lor  the  continuum  equations 
governing  rate-dependent  and  rate-independent  plastic  tlovv  m  metals.  I  he  conser¬ 
vation  laws  are  valid  for  discontinuous  a>  well  as  smooth  solutions  In  the 
rate-dependent  case,  the  evolution  equations  are  in  divergence  lorm.  with  the 
plastic  strain  being  passively  eonveeted  j.'.j  augmented  by  source  terms.  In 
the  rate-independent  case,  the  conservation  laws  involve  a  Lagrange  multiplier 
that  is  determined  by  a  set  of  constraints:  we  show  that  Riemann  problems  tor  this 
system  admit  scale-invariant  solutions.  mo;  Academic  frew  tnc 


1.  Introduction 

In  a  previous  paper  [24],  we  formulated  the  equations  tor  elasticity  in 
the  Euierian  picture  as  conservation  laws.  The  motivation  for  this  work 
was  that  the  Euierian  framework  is  useful  in  numerical  computations  of 
large-deformation  flows.  Furthermore,  the  most  effective  numerical  meth- 
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ods.  such  as  second-order  Godunov  schemes  and  the  trout  tracking  method, 
rely  on  an  understanding  of  the  structure  of  nonlinear  wave  solutions. 
Such  an  understanding  requires  that  the  equations  he  in  conservation 
form.  Subsequent  studies  [5,  I]  have  shown  that  the  use  of  a  conservative 
formulation  does  indeed  lead  to  substantially  improved  numerical  calcula¬ 
tions. 

In  the  present  paper,  we  show  that  a  fully  conservative  formulation  can 
he  given  for  the  equations  governing  an  important  class  of  rate-dependent 
and  rate-independent  plastic  materials.  Unlike  the  case  of  elastic  materi¬ 
als,  for  which  Lagrangian  conservation  laws  were  already  available,  we 
must  appeal  directly  to  experiment  to  decide  which  of  the  dynamical 
equations  for  plastic  How  are  conservation  laws.  For  this  reason,  we  focus 
on  the  plastic  How  of  metals  induced  by  strong  shock  waves.  We  expect, 
however,  that  the  conservative  form  applies  to  a  broader  class  of  plastic 
materials. 

The  formulation  of  Hnite-strain  plasticity  that  we  adopt  has  much  in 
common  with  that  of  Simo  and  coworkers  [30.  28.  29).  There  are  two  novel 
aspects  to  the  present  work.  For  rate-dependent  materials,  we  cast  the 
governing  equations  in  divergence  form  and  argue,  on  the  basis  of  experi¬ 
mental  results  for  metal  plasticity,  that  these  equations  represent  conser¬ 
vation  principles.  For  rate-independent  materials,  we  show  that  the  load¬ 
ing/unloading  conditions,  written  as  an  optimality  condition  involving  a 
Lagrange  multiplier,  hold  for  weak  solutions:  the  resulting  constrained 
system  of  conservation  laws  admits  scale-invariant  solutions. 

Section  2  of  this  paper  contains  a  brief  summary  of  elasticity  in  the 
Lagrangian  and  Eulerian  pictures.  The  various  parts  of  Section  3  present 
the  equations  of  viscoplasticity  in  the  Lagrangian  and  Eulerian  pictures 
and  include  a  discussion  of  the  kinematics  of  plastic  deformation,  constitu¬ 
tive  relations  and  models,  and  plastic  flow  rules.  In  Section  4  it  is  argued, 
on  the  basis  of  experiment,  that  the  plastic  flow  rule  and  hardening  law 
are  conservation  laws.  In  Section  5,  the  conservative  form  for  rate-inde¬ 
pendent  plasticity  is  obtained  from  the  rate-dependent  theory'  in  the  zero 
relaxation-time  limit.  Finally,  in  Section  6,  we  examine  scale-invariant 
solutions  of  the  rate-independent  theory. 


2.  Elasticity 

In  this  section  we  summarize  the  continuum  description  of  an  elastically 
deformable  medium,  first  in  the  Lagrangian  (material)  picture  and  then  in 
the  Eulerian  (spatial)  picture.  Wc  refer  the  reader  to  Ref.  [24]  and  to  the 
general  literature  (e.g.,  Refs.  [15,  35,  21,  9,  2])  for  further  detail. 
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2a.  Kinematics  and  Conservation  Laws 

The  motion  of  a  deformable  body  od  is  represented  mathematically  by  a 
time-dependent  map  </>.  the  motion ,  embedding  $  into  an  ambient  space 
Let  Xn ,  a  =  1.2,3,  denote  material  coordinates  on  od.  and  let  v'. 
/=  1,2.3.  denote  spatial  coordinates  on  then  x‘  -  d>‘{  X.  t ).  The 
motion  cb  must  respect  the  conservation  of  momentum  and  energy'.  These 
principles  are  expressed  as  partial  differential  equations  involving  the 
temporal  and  spatial  derivatives  of  <.6.  namely  the  Lagrangian  velocity 
V1  <>(h‘ /dt  and  the  Lagrangian  deformation  gradient  F[  ■—  <Ub' /i>X" . 

Let  prct  denote  the  mass  density  of  the  undeformed  body,  cf  the  specific 
energy  (energy  per  unit  mass).  Sa,)  the  ( second .  or  symmetric ))  Piola- 
Kirchhoff  stress  tensor,  and  Q"  the  l.agrangian  heat  flux.  Then  the  conser¬ 
vation  iaws  are 


p,X'  =  tc.s'"1)  on 

P.AlKV  -  <f)'-  \V,F‘. $’•),,  -  C>V  0.2) 

Here  the  dot  and  the  semicolon  denote  covariant  differentiation  with 
respect  to  time  and  space,  respectively.  Equations  (2.1)  and  (2.2)  are  in 
divergence  form;  indeed,  they  derive  from  the  integral  form  of  the  conser¬ 
vation  laws.  Therefore  they  hold  in  the  sense  of  distributions. 

Because  the  motion  4>  enters  the  conservation  laws  only  through  its 
derivatives,  the  second-order  equations  can  be  regarded  as  first-order 
equations  in  which  the  velocity  and  the  deformation  gradient  are  funda¬ 
mental  dynamical  variables.  If  <b  is  eliminated  in  favor  of  the  F‘a, 
however,  the  system  of  equations  must  be  expanded.  The  additional 
equations  are  also  conservation  iaws: 

F*  (2-3) 

as  follows  by  equating  mixed  partial  derivatives  of  6.  Equations  (2. 1 )— (2.3) 
express  the  governing  equations  for  a  deformable  medium  as  first-order 
partial  differential  equations  in  conservative  form. 

2b.  Constitutive  Relations 

The  conservation  laws  involve  the  energy,  stress,  and  heat  flux;  there¬ 
fore  they  are  incomplete  unless  supplemented  by  constitutive  relations  for 
these  quantities.  These  relations  characterize  the  material  properties  ot 
the  deformable  medium. 

The  thermodynamic  state  of  an  elastic  material  is  specified  by  the 
deformation  gradient  Fl(t  and  the  specific  entropy  S.  Because  of  covariance 
with  respect  to  rigidbody  transformations  of  the  spatial  coordinates  (the 
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principle  of  frame  indifference),  the  state  depends  on  /•  ',  only  through  the 
(right)  Cauchy -dreen  tensor 


or.  equivalently,  the  Layrunyian  strata  tensor 

- •>.„).  o.5) 

(In  general.  5  (  should  be  interpreted  as  the  material  metric  tensor. )  We 
postulate  that  the  thermodynamic  state  determines  the  specific  energy  6 
through  a  hyperelastic  equation  ol  state 

c  -  ^  (  /i _ V).  (2.6) 


(For  simplicity,  we  assume  that  the  material  is  uniform,  i.e..  neither  /;,cl 
nor  6J  has  explicit  dependence  on  the  material  point  A.)  The  heat  llux  is 
specified  by  a  separate  constitutive  relation,  such  as 

(_)"  =  -.AT",  (2.7) 

where  T  is  the  temperature  and  .A  is  the  thermal  conductivity. 

Following  the  standard  argument  of  Coleman  [3],  the  second  law  of 
thermodynamics,  as  embodied  in  the  Clausius-Duhem  inequality 


iriJ  +  > 

leads  to  the  identifications 


S”*  =  Pro 


M' 


T 


Js  ' 


together  with  the  constraint 


QYTV  >  0. 


2.8) 


(2.9) 

(2.10) 


(2.11) 


For  the  constitutive  relation  (2.7),  this  constraint  requires  that  Jt'>  0.  It 
follows  from  the  governing  equations  and  the  constitutive  relations  that 
the  entropy  satisfies 


Prel™  =  ~Q7:y 


(2.12) 
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Given  the  constitutive  relations  (2.6),  (2.9).  (2.10),  and  (2.7).  Eqs. 
(2.1 )— ( 2.3 )  comprise  a  complete  system  of  evolution  equations  in  conserva¬ 
tive  form  that  describes":elastic  materials  in  the  Lagrangian  picture. 

2c.  Lulerum  Formulation 

To  translate  the  governing  equations  into  conservation  laws  in  the 
Eulerian  picture,  vve  introduce  notation  for  several  Eulerian  quantities, 
which  are  detined  to  he  functions  of  a  -  <f>(  X.  t )  and  t : 


/'  /  '  . 

■J  It  * » 

(2.13) 

(2.14) 

P  :=  1  Wl  • 

(2.13) 

C‘  :=  V. 

(  2.1(7) 

t'  := 

(2.17) 

77  :=  .S’, 

(2.  IS) 

e  ■■=  t. 

(2.19) 

<T,J  ■■=  J-'r^F’e. 

(2.20) 

q‘  :=J-'F\,Q". 

(2.21) 

Here  J  is  the  Jacobian  determinant  of  the  map  <b.  i.e.,  the  determinant  of 
the  linear  transformation  F'a,  and  <r'!  is  the  Cauchy  stress  tensor. 

In  the  standard  fashion,  the  transport  theorem  (see.  e.g..  Ref.  (9])  shows 
that  the  Lagrangian  conservation  laws  (2.1)  and  (2.2)  for  momentum  and 
energy  are  equivalent  to 


jt{pC)  +  {pv V').f  =  rr"  ,,  (2.22) 


-[p(V(r*  4-  e)\  +  [p(kr'  + 


<!k.^ 


2.23) 


The  conservation  law  (2.3)  governing  the  deformation  gradient  trans¬ 
lates  into  the  Eulerian  picture  in  a  similar  manner  [24],  With  the  aid  of  the 
Piola  identity,  {/(T-1)"} ;o  =  0,  Eq.  (2.3)  becomes 

2(0/'.)  +  -  (P'-'/'J:,.  0-24) 


We  supplement  this  continuity  equation  with  the  conservation  of  mass 
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equation. 


<>p 

T7  +  {pr1 ),  -0. 


which  derives  from  die  Lagrangian  equation  pr„,  =  0. 

The  Hulerian  system  of  Eqs.  ( 2.22 )— ( 2.25 )  is  complete  in  the  following 
sense  (ef.  [24]):  Consider  a  solution  to  the  Hulerian  svstem  that  is  either 
smooth  or  piecewise  constant:  suppose  that  there  is  a  deformation  <f>n  such 
that  at  some  initial  time  t  —  f'it  and  ./  are  the  gradient  and  Jacobian  of 
</>„,  respectively.  Then  there  is  a  motion  <l>.  coinciding  with  «/»„  at  t  = 
that  satisfies  the  corresponding  Lagrangian  system  of  Lqs.  ( 2. 1 )— ( 2.3 )  and 
relates  to  the  Hulerian  solution  as  in  the  definitions  above. 

An  alternative  to  Eqs.  (2.24)  and  (2.25)  has  been  suggested  by 
Trangenstein  and  Colella  [34],  viz.,  the  equality  of  mixed  partial  derivatives 
for  <b  ~ 1 : 


0 


=  0. 


(2.26) 


An  equivalence,  in  the  sense  above,  between  the  Lagrangian  system  and 
the  system  comprising  Eqs.  (2.22).  (2.23).  and  (2.26)  is  straightforward  to 
establish  in  the  manner  of  Ref.  [24].  Moreover.  Wagner  [37]  has  used  the 
methods  of  Ref.  [36]  to  extend  the  equivalence  to  general  weak  solutions. 

The  conservation  laws  must  be  supplemented  by  an  equation  of  state, 
which  expresses  the  energy  as  a  function  of  the  Finger  tensor 


6":=/'7\  (2.27) 

or,  equivalently,  the  Almansi  strain  tensor 

;=  4[5;/  -  (6-')„J.  (2.28) 

Note  that  etJ  ■=  (F^lY\En^(F~1)l3r  The  equation  of  state  can  be  written 
[30]  as 


e 


where  r  relates  to  cf  through 


( 2.20) 


*■(««. /\.n)  -  /</>«./*»•  *> 


(2.3(1) 
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Differentiating  this  definition  yields  the  formulae 

--  i>e 

<r‘J  =(>--.  (2.31) 

th’„ 

Oi 

»=—■  (2.32) 

<‘V 

In  calculating  the  derivative  in  Eq.  (2.31).  en  and  fky  are  regarded  as 
independent  variables,  and  their  relationship  is  invoked  only  after  the 
differentiation  has  been  performed. 

Remark.  The  dependence  of  k  on  fky  is  essential.  In  fact. 


P/' 


Ok 

W. 


2(j‘ke 


( 2.33) 


so  that  if  k  were  independent  of  fk y,  then  the  stress  a‘‘  would  vanish.  In 
other  words,  one  cannot  assume  that  the  energy  depends  solely  on  etJ  and 
t]  (unless  the  material  response  to  deformation  is  trivial).  Note,  however, 
that  di/r)fky  vanishes  at  the  undeformed  state  eu  =  0,  which  means  that 
the  effects  of  fky  under  small  deformations  are  of  higher  order.  Similar 
conclusions  hold  in  the  case  of  plasticity. 

The  constitutive  relation  (2.7)  for  the  heat  flux  becomes 


q‘  =  - kd-‘ ,  (2.34) 

with  k  ~  ]~ 

Assuming  the  constitutive  relations  (2.29),  (2.31),  (2.32),  and  (2.34),  Fqs. 
(2.22)  and  (2.23).  together  with  either  Eqs.  (2.24)  and  (2.25)  or  Eq.  (2.26). 
comprise  a  complete  set  of  conservation  laws  in  the  Eulerian  picture. 


3.  Plasticity 

The  plastic  response  of  real  materials  is  extremely  diverse.  For  instance, 
the  microphysical  pictures  of  metal  plasticity  and  of  the  flow  of  granular 
materials  are  quite  different.  Furthermore,  passage  from  the  microphysical 
to  the  macroscopic  description  requires  scientific  modeling  across  a  hierar¬ 
chy  of  length  scales.  Nevertheless,  there  are  two  concepts  that  apply  to  a 
broad  class  of  macroscopic  models:  plastic  deformation  and  plastic  flow 
rule. 
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3a.  Plastic  Deformation 

One  principal  feature  p.f  plastic  behavior  is  irrecoverable  deformation: 
whereas  an  elastic  body  will  return  to  its  undeformed  state  when  surface 
forces  are  relaxed,  a  plastic  material  might  not.  The  permanent  deforma¬ 
tion  suffered  depends  on  the  manner  in  which  the  body  was  worked,  i.e.. 
on  the  history'  of  deformation  of  the  material.  In  a  continuum  theory'  of 
this  phenomenon,  there  must  be  “internal'’  state  variables  that  record  this 
deformation  history.  Because  the  deformation  history  varies  from  point  to 
point  in  the  body,  a  plastic  material  is  inhomogeneous. 

To  understand  these  ideas  on  a  physical  level,  consider  a  small  neigh¬ 
borhood  of  a  material  point  T  in  a  body  that  has  undergone  plastic 
deformation.  If  the  neighborhood  is  cut  out  of  the  material  and  the  forces 
on  its  surface  are  relaxed,  then  it  deforms,  recovering  a  stress-free  config¬ 
uration.  Because  of  the  plastic  working,  this  configuration  differs  from  the 
original  configuration  of  the  neighborhood  before  the  body  was  deformed. 
Thus  at  every  point  there  is  a  local  irrecoverable  deformation  taking  the 
neighborhood  from  its  virgin  state  to  a  stress-free  configuration.  In  gen¬ 
eral,  however,  there  is  no  global  deformation  that  relaxes  the  internal 
stresses  simultaneously  throughout  the  body.  The  local  irrecoverable  de¬ 
formations  at  neighboring  points  may  be  incompatible,  in  the  sense  that 
the  relaxed  neighborhoods  do  not  fit  together  to  form  a  continuous  body 
without  gaps  or  tears. 

These  ideas  can  be  expressed  mathematically,  using  the  theory  of 
inhomogeneous  materials  developed  by  Noll  [35].  For  a  neighborhood  that 
is  sufficiently  small  on  the  macroscopic  scale,  the  irrecoverable  deforma¬ 
tion  is  approximated  well  by  its  gradient,  which  is  a  linear  map  from  the 
tangent  space  at  X  to  three-dimensional  Euclidean  space  R\  Mathemati¬ 
cally,  such  a  linear  map  defines  a  frame  at  X,  for  its  inverse  carries  the 
standard  basis  for  R3  to  a  basis  for  the  tangent  space  at  X.  We  denote  the 
components  of  this  local  reference  frame  at  X  by  ( Frcf(  X  )Ya ,  where 
a  =  1,2,3. 

Let  a  local  reference  frame  (Fref(30)"a  be  chosen  at  each  material  point 
X.  As  discussed  above,  there  need  not  exist  a  deformation  of  the  entire 
body  whose  gradient  is  (Fref)aa:  there  is  no  “intermediate”  stress-free 
configuration  of  an  inhomogeneous  body.  Indeed,  the  body  is  considered 
to  be  homogeneous  if  and  only  if  the  frame  field  Fra  is  a  gradient. 
Nevertheless,  the  gradient  F‘a  of  an  arbitrary  deformation  can  be  re¬ 
ferred,  at  each  material  point,  to  the  local  reference  frame.  Thus  we 
define  F'a  by  decomposing  the  deformation  gradient  F\t  as  a  product 
involving  the  frame  (FTef)aa: 


F‘a(Fra) 


a 

a  * 


(3.1) 
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Fio.  1.  The  multiplicative  decomposition  of  'he  deformation  gradient:  F‘n  =  ( Fe YJ  FpY‘n. 
The  plastic  deformation  (Fp)“„.  which  represents  the  irrecoverable  deformation,  maps 
neighborhoods  in  the  body  to  local  stress-free  configurations.  Here  this  deformation  involves 
the  slip  of  horizontal  crystal  planes.  As  indicated,  the  relaxed  neighborhoods  do  not  ht 
together  continuously.  The  elastic  deformation  (Fc)'u  distorts  the  relaxed  neighborhoods  into 
the  final  configuration  of  the  body. 


(See  Fig.  1.)  For  a  homogeneous  body,  the  transformation  to  local  refer¬ 
ence  frames  amounts  to  a  change  of  material  coordinates. 

For  example,  the  frame  field  can  be  given  a  concrete  interpretation  in 
terms  of  the  theory  of  dislocations  in  crystals  (see.  e.g.,  Refs.  [35,  6,  13.  12. 
10]).  Each  of  the  three  basis  vectors  defined  by  the  local  reference  frame 
at  X  is  regarded  as  a  displacement  by  a  fixed  number  of  lattice  spacings 
along  a  crystal  axis.  Thus  the  local  reference  frames  reflect,  at  a  contin¬ 
uum  level,  the  crystallographic  structure.  The  tensor  field 

P”;=  -£"”•(  Fref)%:,,  (3.2) 

which  is  a  measure  of  the  nonintegrability  of  the  frame  field,  can  be 
interpreted  as  the  dislocation  density. 

With  ( Frel chosen,  the  material  response  depends  on  F(J,  the  en¬ 
tropy,  and  the  material  position  X.  If  there  exists  a  choice  of  local 
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reference  frames  sueli  that  this  response  does  not  depend  explicitly  on  X. 
then  the  body  is  said  to  he  materially  uniform.  Note  that  there  is  a  possible 
ambiguity  in  the  decomposition  (3.1).  depending  on  the  symmetry  of  the 
material,  indeed,  a  material  isomorphism  at  a  point  is  a  change  of  local 
reference  frame  that  preserves  the  response  of  the  material  to  deforma¬ 
tions  [35.  K)].  These  isomorphisms  form  the  local  symmetry  group  at  the 
point:  a  materially  uniform  body  has  the  same  local  symmetry  group  at 
each  point.  One  is  tree  to  transform  the  local  reference  frame  at  each 
point  of  the  body,  provided  that  the  transformation  belongs  to  the  symme¬ 
try  group  and  depends  continuously  on  position.  Fur  a  crystal,  the  symme¬ 
try'  group  is  finite,  so  that  the  transformation  must  Ire  the  same  at  each 
point.  Thus  ( l'lctY\,  is  essentially  unique,  being  defined  modulo  a  finite 
group  of  transformations.  By  contrast,  the  symmetry  group  is  the  rotation 
group  lor  a  material  such  as  a  metal.  Consequently,  the  constitutive 
equation  tor  the  energy  depends  on  only  through  the  strain 


(3.3) 


The  theory  of  inhomogeneous  materials  can  be  adapted  to  model 
plasticity  by  allowing  the  local  reference  frames  of  Noll  to  be  dynamic  [6], 
In  this  approach,  the  frame  field  ( /’Ic,  )'„  is  regarded  as  the  plastic 
deformation  (Fp)"(t.  while  is  the  elastic  deformation  (/•[.)',.  Thus  we  are 
led  to  the  fundamental  kinematic  decomposition 


r..  =  (Fv)'jrt 


(3.4) 


Such  a  multiplicative  decomposition  has  been  advocated  by  several  work¬ 
ers  in  plasticity  theory  [4,  17.  6.  16.  41.  20.  16.  25.  14j. 

The  plastic  deformation  ( Fp )"rt  is  the  principal,  but  not  necessarily  the 
only,  internal  variable  in  a  model  of  a  plastic  material.  Other  variables, 
such  as  a  hardening  parameter ,  characterize  the  effects  of  plastic  deforma¬ 
tion  on  the  material  response.  We  will  assume  that  there  is  a  single  scalar 
hardening  parameter  K\  the  theory  is  easily  extended  to  account  for 
several  tensoriai  internal  variables. 


3b.  Constitutive  Relations 

Internal  variables  (Fp)"„  and  K  enter  the  equations  of  motion  in  two 
places:  in  the  dynamical  equations  that  govern  their  evolution;  and  in  the 
constitutive  equation  for  the  energy,  which  depends  on  the  local  reference 
frame  and  hardening  parameter  as  well  as  the  deformation  gradient  ami 
the  entropy.  Choosing  a  plasticity  model  amounts  to  specifying  the  evolu¬ 
tion  and  constitutive  equations. 
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A  general  thermodynamical  framework  for  plasticity,  which  encom¬ 
passes  a  variety  of  specific  models,  has  been  given  by  Green  and  Naghdi 
[7].  /his  framework  allows  for  various  choices  of  the  measure  of  plastic 
strain.  A  model  based  on  the  multiplicative  decomposition  (3.4)  can  be 
accommodated  in  two  ways,  depending  on  the  material  symmetry'  [6],  For 
crystals,  one  chooses  (/rp)‘,ii  as  the  plastic  strain  measure;  for  metals,  one 
defines  the  plastic  strain  Fvtfi  by  Eq.  (3.3); 

c:;,  -  o-n 

in  the  following  we  will  focus  on  metals,  although  the  development  could 
be  adapted  to  treat  crystals. 

The  constitutive  equation  for  the  energy  expresses  ef  as  a  function  of 
the  deformation  gradient  F't,  the  plastic  deformation  (F  )"„,  the  harden 
ing  parameter  K.  and  the  entropy  S.  By  the  principle  of  frame  indiffer¬ 
ence,  £  depends  on  the  deformation  gradient  only  through  the  (total) 
strain  Enji\  and  for  a  metal,  the  plastic  deformation  appears  only  through 
the  plastic  strain  Therefore  we  assume  the  hyperelastic  equation  of 
state 

-?(£*,  E’s,K.S).  (3.6) 

This  assumption,  together  with  the  Clausius-Duhem  inequality,  implies 
formulae  for  the  stress  and  temperature  [3,  7]. 

To  derive  these  formulae,  we  introduce  the  notation 

+  (3.7) 

IK,  “  -  P«;.F*»)  ( 3-8) 

for  the  Lagrangian  strain  rate  and  rorticity  tensors.  These  tensors  are  the 
symmetric  and  antisymmetric  parts  of  Fl.aFkf),  respectively;  in  particular, 

F„ti  ~  Dap-  (3.9) 

Note  that  the  energy  conservation  equation  (2.2),  con.bined  with  the 

momentum  conservation  equation  (2.1),  implies  that 

-  <2V  (3.10) 

Furthermore,  differentiating  the  equation  of  state  (3.6)  shows  that 
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Combining  these  two  equations  with  inequality  (2.8)  yields 


5  P  /Cel  ti  /'re  I  -  p  /)rct  ^ 


/V, 


QyT, 


3.12) 


Assuming,  as  is  conventional,  that  Etl0  and  S  can  be  chosen  arbitrarily 
even  if  there  is  neither  a  temperature  gradient  nor  plastic  flow  (i.e.. 
T  =  0.  -  U.  and  K  =  0).  inequality  (3.12)  leads  to  the  identifications 


=  o 

•’  /'ret 


(3.13) 


along  with  the  constraints 


7  =  — 
r/5 


27  >  0. 


-  +  i?  >  o. 


(3.14) 


(3.15) 

(3.16) 


Here 


^  Pre(dE^E"ti  Prc>flKK 


(3.17) 


is  the  plastic  dissipation .  since  entropy  is  generated  according  to 

Pref75-  -<2*:y  + 

Thus  the  piuStic  dissipation  is  3  =  +  UK.  where 


(3.18) 


‘/Act  ,r-p 


(3.19) 


n  :  Pref, m 


(3.20) 


are  the  dissipative  forces  conjugate  to  E*p  and  K.  respectively.  In  general. 
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the  plastic  stress  .SI’*1  end  the  stress  .S' "l>  ditier,  although  they  are  approxi¬ 
mately  equal  in  the  small-strain  limit.  Inequalities  (3.15)  and  (3.16)  place 
thermodynamic  restrictions  on  the  dependence  ot  6  on  the  internal 
variables  and  on  the  evolution  equations  that  govern  them. 

3c.  Constitutin'  Models 

Specifying  the  constitutive  funen  >n  6  involves  modeling  assumptions. 
For  instance,  a  natural  assumption  in  plasticity  is  that  the  energy  is 
determined  primarily  by  the  elastic  part  of  the  strain:  the  history  ot  plastic 
deformation  has  only  mild  effects  on  the  material  response,  accounted  for 
by  the  hardening  parameter,  furthermore,  the  anisotropic  component  ot 
the  elastic  strain  is  typically  small  because  ot  plastic  yielding,  so  that  a 
small-anisotropy  approximation  is  useful  for  a  material  that  is  isotropic  in 
the  local  stress-free  configurations. 

Two  alternative  measures  of  the  elastic  strain  are  used  in  the  literature. 
Green  and  Naghdi  define  the  elastic  strain  as  measured  with  respect 
to  the  Lagrangian  frame,  through  the  additive  decomposition 


of  the  Lagrangian  strain.  Equation  (3.21)  is  analogous  to  the  additive 
decomposition  of  infinitesimal  strain  in  classical  plasticity.  Alternatively. 
Lee  and  Liu  [17.  16]  measure  the  elastic  strain  with  respect  to  the  local 
reference  frames  using 

C  -  -  «,„]•  0-221 

instead  of  E^.  These  strain  measures  are  related  by 

El?  =  (3.23) 

Either  measure  of  elastic  strain  can  be  used  in  constructing  the  equation 
of  state  (3.6).  This  is  evident_for  the  additive  decomposition  (3.21).  In 
terms  of  the  strain  measure  E^,  the  general  constitutive  equation  (3.6) 
can  be  written  as 

(3.24) 

where 

*(&.(*;)*„.  K.s)  )“„£,%( /-p )", ,  +  i[(F„ . 

■- 1  ('vU'vL, |.K..q. 


(3.25) 
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Conversely,  for  Eq.  (3.24)  to  represent  an  equation  of  state.  Eq.  (3.25) 
must  hold  for  some  £ ’.  Indeed,  when  the  local  symmetry  group  is  the 
rotation  group.  E*h  and  ( fp)"„  enter  6’  only  through  the  variables  Ell(i 
and  Eptli.  In  this  sense,  the  alternative  formulations  are  equivalent  [8,  27]; 
for  a  particularly  lucid  explanation,  see  Ref.  [30]. 

When  the  equation  of  state  is  formulated  as  in  Eq.  (3.24).  it  is  useful  to 
work  with  various  tensors  referred  to  the  intermediate  frame.  For  in¬ 
stance.  with 


•?**-  O'p (3.26) 
Eqs.  (3.13)  and  (3.25)  imply  that 

Sab  =  Fro  ~=^r  ■  (3.27) 

Furthermore,  in  formulating  the  dynamical  equation  governing  the  plastic 
strain  (see  Section  3d),  the  following  definitions  will  be  helpful; 

(3.28) 

{L’i'U-  (3.29) 

Here  ( LpEv)ilh  is  the  plastic  Lie  derivative  of  Eph  [30];  cf.  Eq.  (3.63)  below. 

Let  us  return  to  the  example  of  a  plastic  material  that  is  isotropic.  The 
physical  picture  provided  by  the  multiplicative  decomposition  of  the  defor¬ 
mation  gradient  (see  Fig.  1)  suggests  that  the  isotropy  is  manifested  in  the 
stress-free  intermediate  configuration.  This  means  that  <T_depends  solely 
on  the  principal  invariants  of  (EcYh,  or  equivalently  of  (Cc)ah  =  2 (Ee)ah 
4-  8“h.  (Recall  that  the  principal  invariants  of  a  3x3  matrix  A  are 
i}(A)  =  trA.  i2(A)=\[(lrA)1-{rA1],  and  t,(xf)  =  det  A:  see,  e.g.. 
Refs.  (2,  9].)  Thus  Lee  and  Liu  [17.  16]  adopt  the  equation  of  state 

<f  =  <f(  AT,  S),  (3.30) 

where  Ik  ■—  iA.(Ce),  k  =  1,2,3.  This  equation  of  state  conforms  with  the 
requirement  embodied  in  Eq.  (3.25)  because  the  invariants  /,.  /-,.  and  /, 
can  be  expressed  in  terms  of  and  Ep(j.  To  see  this,  we  let  = 
2Epp  +  8af)  and  note  that 


(Q'f  -  K-T.tc')  V,)V 


(3.31) 
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as  thejnvariants  are  unaffected  by  a  similarity  transformation,  ii(C~,C) 
=  it(CcU  =  l,2,3. 

Remark.  For  an  equation  of  state  such  as  Eq.  (3.30),  the  plastic  stress 
5"^  of  Eq.  (3.19)  that  enters  the  plastic  dissipation  _?  can  be  related  to  the 
stress  S"li: 


(the  matrix  on  the  right  is,  in  fact,  symmetric).  When  the  elastic  strain 
is  small,  is  approximately  equal  to 

3d.  Plastic  Plow  Rule 

The  other  fundamental  feature  of  plastic  behavior  is  the  plastic  flow 
rule ,  which  specifics  the  time  evolution  of  the  plastic  deformation.  Such  a 
How  field  has  the  general  form  [23.  18] 

=  \ali(Eyti.E*A.K.S)  (3.33) 

for  a  rate-dependent  plastic,  i.e..  riscoplastic.  material.  In  addition,  a 
hardening  law 

K  =  H(Ey/i,  ££>,  K.S)  (3.34) 

governs  the  evolution  of  the  hardening  parameter.  The  flow  rule  and 
hardening  law  are  constrained  by  the  requirement  (3.15)  that  the  plastic 
dissipation  3  be  nonnegative. 

Characteristic  of  viscopiastic  flow  (as  opposed  to  viscoelastic  flow)  is 
that  plastic  deformation  occurs  only  if  a  threshold  has  been  reached. 
According  to  the  plastic  flow  rule  (3.33),  the  material  deforms  plastically 
only  in  the  region  of  state  space  where  Aafi  #  0.  Since  the  hardening 
parameter  records  the  effects  of  plastic  deformation  on  the  material 
properties,  hardening  occurs  only  during  plastic  flow;  therefore  H  is 
assumed  to  vanish  when  Anfj  =  0.  The  interior  of  the  region  where 
Aafj  -  0,  in  which  the  behavior  is  purely  elastic,  is  called  the  elastic  range. 
The  boundary  of  the  elastic  range  is  called  the  ( static )  yield  surface.  The 
elastic  range  and  yield  surface  are  usually  characterized  by  a  yield  function 
<t>:  <  0  in  the  elastic  range  and  4>  =  0  on  the  yield  surface.  The  yield 

function  <I>  depends  on  the  variables  E.s,  Efs,  K.  and  S:  frequently, 
however,  <t>  depends  on  the  strains  and  the  entropy  only  through  the  stress 
and  the  temperature. 

We  illustrate  these  concepts  with  a  model  that  has  been  used  to 
describe  high  strain-rate  plastic  flow  in  metals  [11,  31,  32],  In  the  case  of 
finite  deformations,  this  model  is  specified  most  simply  in  terms  of  tensors 
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referred  to  the  intermediate  frame  (see  Section  3c  and  the  first  remark 
below).  First,  the  ran  Mises  yield  criterion 

<!>  :=  lidev  .VII  -  /273  T„(  F.  T.  K  )  ( 3.35) 

is  adopted.  In  this  equation.  Y()  is  the  static  yield  strength  (as  specified 
below).  P  ■-  \J~XS\  is  the  mean  pressure,  (dev  S )“h  is  the  deriatoric 
stress. 


( dev  S)'‘h  =  S‘,h  -  \S‘\.  d“h.  ( 3.36) 

and  the  norm  jj  •  II  is  defined  by  11.411'  =  A“hA  Jh.  The  flow  ruie  is  taken  to 
be  the  Lecy-St.  Venant  flow  ride,  given  by 


-  (dev  S)llh 

(Epfp),„,=  \(r.WdcvS\\,T,K)  (3.37) 

Ildev.SII 

In  particular,  the  model  is  complete  once  ||LP£P||  =  \(P.  ||dev  5||.  T.  K ) 
has  been  specified. 

Remarks.  (1)  Plastic  flow  in  metals  is  often  assumed  to  be  volume 
preserving  [17.  16],  so  that  det  Fp  =  1.  Since 

(dec  F„)  -  '(del  f„)  -  (  Fr)\(  F~  (3.38) 

the  plastic  flow  ruie  preserves  volume  if  and  only  if  (LpE)ah  is  trace-free; 
this  has  been  emphasized  to  us  by  F.  Wang  [40],  For  this  reason,  we  have 
formulated  Eqs.  (3.35)  and  (3.37)  in  terms  of  (dev.V)ol\  rather  than 
(devo-)",  as  would  be  conventional  (see,  e.g.,  Ref.  [28]).  The  correspond¬ 
ing  Eulerian  formulation  is  given  in  Eqs.  (3.71)  and  (3.72)  below. 

(2)  When  Y{)  is  independent  of  F ,  the  Levy-St.  Venant  flow  rule  is 
associative ,  in  the  sense  that 

r?<J> 

{  LpEp).h  -  A(/>,lldev5ll,r,K)~=^.  (3.39) 

Thus  the  direction  of  \a(i  is  determined  by  the  normal  to  the  yield  surface 
in  stress  space. 

(3)  In  terms  of  Lagrangian  quantities,  the  flow  rule  (3.37)  is 


A(  P,  ||devrp  Sllc-p,  T,  AT)Cpy 


(devfP.V)yd 

||devrP5||rpC^ 


(3.40) 
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where  P  =  \J~\SyiC^, 

(d  cv(;>  S  )'*■■=  S,,(i  -  \SynC^(Cp  (3.4!) 

and  the  norm  j|  •  ||(t  is  defined  hv 


Therefore  t his  flow  rule  does,  indeed,  have  the  form  of  Lq.  (3.33). 
provided  that  .S'"p  and  /'  are  expressed  in  terms  of  K .  and  S 

through  Eqs.  (3.13)  and  (3.14). 

The  specification  of  ||/J7:P!|  derives  from  a  model  for  the  dynamic  yield 
strength  Y  —  T(  P.  T.  ||/J7-.'P1I.  K ).  1'his  quantity  is  the  yield  strength  ob¬ 
served  in  plastic  flow  at  nonzero  strain  rate:  in  other  words. 

Ildev  .VII  =  VV3  Y(  P.  T.  \\LpEp\\.  K  )  (  3.43) 

during  plastic  flow.  The  dynamic  yield  strength  must  exceed  the  static  yield 
strength 


Ya(  P.  T .  K  )  -  Y(  P,  /',  0.  K  )  ( 3.44) 


and  increase  monotonically  with  strain  rate. 

Given  a  model  for  the  dynamic  yield  strength,  the  corresponding 
flow  rule  is  derived  as  follows.  When  Ildev  .V l|  is  below  the  static  yield 
strength,  there  is  no  plastic  flow,  and  IILP/:P||  =  0:  otherwise.  Ildev  ,V||  > 
y/2/3  Yn(  P.  T.  K).  so  that  f£q.  (3.43)  can  be  solved  for  ||/.p/Tp|i  in  terms  of 
the  remaining  quantities,  obtaining 

||LpTp|l  =  A (/\  Ildev  .VII.  7\  K).  (3.45) 

Therefore  the  plastic  flow  rule  is 


0  if  Ildev  .VII  <  v 2/3  V()(  P.  T.  K). 

-  -  -  (dev .V)„/5 

A(  P,  ||dev  .V||,  T,  K ) - rr—  otherwise. 


To  be  more  concrete,  we  consider  a  model  for  the  dynamic  yield 
strength  given  by  Steinberg  and  Lund  [32],  In  this  model,  Y  is  composed  of 
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a  thermally  activated  part  )',  and  a  strain-hardening  part  )',: 

(i(P.  /’), 

V  ^  - !  v,  -  V,(  K)  ].  (3.471 


where  the  thermally  activated  part  Y, .  0  <  V,  <  Yr.  is  the  solution  o!  the 
equation 


/."/•  |,!  -  Y . 


(  . 


-  (  ,  ' Y ,  C.\p 


Y, 

Yr 


( 3.4<S) 


In  these  expressions,  the  lunctions  (l  and  V,  and  the  parameters  (,. 
C\.  ))..  and  l  \/k  are  mtitenal-dependent  quantities.  The  hardening  law 
(3.34)  is  taken  to  be 

K  =  v  2/F[\L,'Frl\:  (3.49) 

in  other  words.  K  is  the  equivalent  plastic  strain,  and  Eq.  (3.49)  models 
isotropic  strain-hardening . 

The  function  A  for  this  model  is  obtained  as  follows:  According  to  Eq. 
(3.48),  Y,  =  d  when  i|/J7:r:!  =  0.  so  that 

G(T\T) 

>’„(  r.  r.  k  )  =  — - —  y,(  K).  ( 3.50) 

*7n 

Since  ||dev.V‘i  =  ^2/3  Y  during  plastic  How.  Eq.  (3.47)  implies  that 

I  dev  .VII  -  yjl/F  V„(  F.  7  .  K  ) 

y  =  ___ —  - =- - - :  (3.51) 

V'2/3G(/E  /  )/(;„ 

thus  Y,  coincides  with  a  measure  of  the  extent  to  which  the  static  yield 
condition  has  been  exceeded.  Substituting  this  expression  into  the  right- 
hand  side  of  Eq.  (3.48)  gives  ^2/ 3  A. 


3c.  Lulerian  Formulation 

Corresponding  to  the  internal  variables  and  K  characterizing  a 
plastic  material  in  the  Lagrangian  picture,  we  define  the  Eulerian  quanti¬ 
ties 


/>  i1 


(T-1) 


K  :=  K, 
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which  are  regarded  as  functions  of  .r  =  d>(  X.  i )  and  /.  In  terms  of  the 
elastic  Finger  tensor 


K' 

(3.55) 

(3.56) 

The  equation  of  state  is 

t-  * 

(3.57) 

where 

*(»«.«*. /%.«.*)  -  <H/'V . r,./VS./*4. *•*?).  (Ml 


Just  as  for  elastic  materials,  the  Eulerian  conservation  laws  of  momen¬ 
tum.  energy,  and  continuity  are  Eqs.  (2.22).  (2.23),  and  either  Eqs.  (2.24) 
and  (2.25)  or  Eq.  (2.26);  the  constitutive  equations  are  Eqs.  (3.57).  (2.31). 
(2.32),  and  (2.34).  The  thermodynamic  constraints  (3.15)  and  (3.16)  be¬ 
come 


d  >  0.  (3.59) 

qk6  k 

- ~-  +  d  >  0.  (3.60) 

where 

de  r>£ 

d  :=  -p  —  (L,ep)tl  -  P  —  L,k  (3.61) 

de?:  ok 


is  the  plastic  dissipation  J~lS2,  and  entropy  is  generated  according  to  the 
equation 

pdLtrf  —  —qkk  +  d.  (3.62) 

In  these  equations  we  have  introduced  the  Lie  derivatives  [21,  30]  of  e?r  k. 
and  t],  as  defined  by 

(L,  e-),j  - 


=  +  vk*U  +  vk:A,  +  efck:r 

(3.63) 

Sk 

L,.K  ;  = - h  VkK.  , 

'  dt 

(3.64) 

dr> 

L..1J  := - 1-  V  7i.k. 

L 1  dt  ,k 

(3.65) 
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To  write  the  plastic  How  rule  (3.33)  in  the  Eulerian  formulation,  we 
multiply  Eq.  (3.33)  by  ( /' 1  >"(  /  1  )0,  to  find  that 

(L.ep)n  =  A  (3.66) 

in  addition,  the  hardening  law  (3.34)  becomes 

L:K  =  fl{eu.e'kl,  fky,K.ri).  (3.67) 

Here  we  have  made  the  identifications 


'V,  {L'kl'  Ckl'  J  ' 


.*■■?)  =  ( c 1  r,A..»(  tv-  -  *  )<'•"’  >',• 


(3.68) 

(3.69) 


Similarly,  the  Eulerian  yield  function  is  given  by 

*(*«•<*./*,.«•»»)  =<!>(£•„.  E!,.K.  S).  (3.70) 

For  example,  the  von  Mises  yield  criterion  (3.35)  and  the  evy- 
St.  Venant  flow  rule  (3.37),  as  expressed  in  the  Eulerian  picture,  are 


3.71) 


,h  (devour) 


ki 


(3.72) 


< p  l|dev,,;  i<rit/v  i  -  y'2/3  y.,(  p.  0,  k). 

Lre,p  =  A(  /;,  ildev,,  - 1  <r||/,  ■  ,  0,  k)(6c~  1  )ik 

1  v  c  1  7  ||devh-  ltrii^,  < 

where  p  ■=  \(rkl(b~  1  )kl  is  the  mean  pressure  P. 

(devour)"  -<r"  -  (3-73) 

and  the  norm  ||  •  IU-i  is  defined  by 

\\A\\l^-=A‘>{bL:')lkAk,{b^),r  (3-74) 


Also,  A  and  v„  are  identified  with  A  and  y(),  respectively.  Therefore  this 
flow  rule  has  the  form  of  Eq.  (3.66)  when  a'*  and  0  are  expressed  in  terms 
of  ekl,  e%„  k,  and  77  through  Eqs.  (2.31)  and  (2.32). 


4.  Flow  Rules  as  Conservation  Laws 

In  the  previous  section  we  have  seen  that  the  evolution  of  a  plastic 
medium  is  governed  by  the  conservation  laws  of  momentum,  energy,  and 
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continuity  together  with  flow  rules  for  the  internal  variables.  This  discus¬ 
sion  implicitly  assumed  that  the  flow  was  smooth.  Here  we  consider  the 
possibility  of  discontinuous  solutions,  as  are  expected  to  occur  on  the  basis 
of  experiment  and  the  structure  of  the  governing  system  of  partial  differ¬ 
ential  equations.  Analysis  of  discontinuous  solutions  is  best  carried  out 
when  the  equations  are  written  as  conservation  laws.  This  poses  a  ques¬ 
tion:  When  are  the  dynamical  equations  for  the  internal  variables  conser¬ 
vation  laws? 

Answering  the  question  of  when  an  evolution  equation  is  a  conservation 
law  has  two  aspects.  First,  the  equation  must  be  written  in  divergence 
form.  Second,  the  equation  must  hold  in  the  sense  of  distributions.  To 
understand  these  aspects,  we  recall  the  example  of  gas  dynamics. 

For  smooth  gas  flows,  the  equations  governing  the  time  evolution  of 
both  energy  e  and  entropy  17  are  in  divergence  form: 

-[p(ir,r'  +  e )]  +  [pG'T'  +  =  ~{prJ):n  (4.1) 

d 

—  (pT?)  +  (P*7'')w  =  0.  (4.2) 

However,  one  cannot  require  both  of  these  equations  to  hold  for  discon¬ 
tinuous  flows.  Indeed,  the  jump  conditions  enforcing  conservation  of 
energy  entail  an  increase  in  entropy,  whereas  the  jump  conditions  for 
conservation  of  entropy  lead  to  a  decrease  in  energy'.  Although  both  Eq. 
(4.1)  and  Eq.  (4.2)  hold  up  until  shock  formation,  a  choice  between  them 
must  be  made  at  that  point.  This  choice  is  not  dictated  by  the  mathemati¬ 
cal  structure  of  the  equations.  Rather,  it  must  be  made  on  the  basis  of 
physics:  energy  is  conserved  and  entropy  increases.  Similarly,  in  plasticity, 
we  must  appeal  to  physics  to  decide  when  a  flow  rule  is  expressing  a 
conservation  law. 

For  plasticity,  we  first  ask  whether  the  evolution  equations  for  the 
internal  variables  have  divergence  form.  Referring  to  Eqs.  (3.33)  and 
(3.34),  we  see  that  the  plastic  flow  rule  and  hardening  law  have  divergence 
form  in  the  Lagrangian  formulation:  the  time  derivatives  of  the  internal 
variables  are  given  by  source  terms  involving  the  state  variables  but  not 
their  derivatives.  By  contrast,  the  corresponding  equations  (3.66)  and 
(3.67)  in  the  Eulerian  formulation  have  not  been  written  in  divergence 
form;  this  is  because  L,.ep  and  Lt  k  involve  spatial  gradients  multiplied  by 
functions  of  the  state  variables.  We  therefore  replace  these  equations  by 
ones  in  divergence  form. 

The  Eulerian  equations  in  divergence  form  that  correspond  to  the 
Lagrangian  flow  rule  and  hardening  law  are  obtained  as  in  Ref.  [24]: 
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(a)  the  equations  are  written  in  integral  form:  (b)  a  change  of  variable  is 
made  in  the  integral;  and  (c)  the  transport  theorem  is  applied.  The 
equations  that  result  can  also  be  obtained  by  a  simple  argument  (which  is 
valid  only  for  smooth  solutions).  Thus  we  multiply  Eq.  (3.33)  by  p  and  add 
the  equation  of  conservation  of  mass,  multiplied  by  this  yields 


iH 


(p*U)  * 


(4.3! 


Since  'vc  obtain 


+  ( w..  <■,",/ ).» -  watv  (4.4) 


Similarly,  the  hardening  law  becomes 


—  {px)  +  (pKrk).k  =  ph.  (4.5) 

at 

The  question  remains  whether  Eqs.  (4.4)  and  (4.5)  represent  physical 
conservation  laws.  In  this  paper,  we  address  this  question  within  the 
context  defined  by  experiments  on  shock  wave  propagation  in  metals  [26. 
38,  39,  31-331. 

Distinct  wave  structures  are  observed  in  the  three  regimes  of  low. 
intermediate,  and  high  pressure  driving  the  metal.  If  the  driving  pressure 
is  below  the  Hugoniot  elastic  limit  (on  the  order  of  10  kbar).  a  purely 
elastic  shock  wave  occurs.  For  an  intermediate  pressure  (less  than  about 
100  kbar),  a  two-wave  structure  is  observed.  At  still  higher  pressures,  the 
two-wave  structure  collapses  into  a  single  stron^  shock  wave. 

In  the  intermediate  pressure  regime,  the  leading  wave  is  relatively 
sharp,  has  an  amplitude  that  is  generally  independent  of  the  driving 
pressure,  and  propagates  at  a  speed  determined  by  the  elastic  properties 
of  the  metal.  Therefore  it  is  interpreted  as  an  elastic  precursor  wave.  The 
wave  following  the  precursor  has  a  broader  profile  and  an  amplitude  that 
increases  with  driving  pressure;  its  structure  is  sensitive  to  the  strain  rate. 
Furthermore,  the  relaxation  time  associated  with  the  strain-rate  depen¬ 
dence  is  long  compared  with  the  rise  time  for  the  elastic  precursor.  These 
observations  suggest  that  the  plastic  deformation  occurs  predominantly 
after  the  passage  of  the  elastic  precursor. 

On  this  basis,  it  is  concluded  that  (a)  the  elastic  wave  can  be  modeled 
effectively  as  a  mathematical  discontinuity,  involving  no  change  in  the 
plastic  deformation  and  (b)  the  plastic  flow  occurs  within  a  broad  profile 
plastic  wave  following  the  elastic  precursor.  This  picture  is  consistent  with 
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the  assumption  that  the  flow  rule  and  hardening  iaw.  Eqs.  (3.33)  and  (3.34) 
or  Eqs.  (4.4)  and  (4.5),  are  genuine  conservation  laws. 

To  justify  this  statement,  we  consider  a  plane  wave,  propagating  at 
(Lagrangian)  speed  Sv  in  the  direction  of  the  vector  Nit,  across  which  a 
discontinuity  occurs  in  the  state  variables.  Using  A  to  denote  the  jump  in  a 
quantity  across  the  wave,  the  jump  conditions  for  Eqs.  (2.1 )— (2.3)  and  Eqs. 
(3.33)  and  (3.34)  are 


=  A(F,,S‘*).V  (4.6) 

-/>rc,AvA(iKr  +  <f)  =  A(I/F,.V^)%  -  ±Q>Ny.  (4.7) 

-5.VAF'„  =  A  V'Nt.  (4.8) 

-  A v  A E%p  =  0.  (43)) 

-SNAK  =  i).  (4.10) 


Therefore  no  change  in  the  internal  variables  Ep  and  K  are  sustained 
unless  SN  =  0;  i.e.,  the  wave  moves  with  the  speed  of  the  fluid  particles. 

Thus  for  the  case  of  shock  wave  propagation  in  metals,  we  have 
established  a  conservative  form  for  the  governing  equations.  This  conser¬ 
vative  form  might  well  apply  to  a  broader  class  of  rate-dependent  plastic 
materials. 


5.  Rate-Independent  Plasticity 


While  plastic  flow  is  fundamentally  a  rate-dependent  process,  a  rate- 
independent  approximation  is  useful  in  certain  flow  regimes.  In  this 
section,  we  shall  see  how  the  rate-independent  equations  emerge  as  a 
singular  limit  of  rate-dependent  plasticity.  We  also  explain  the  sense  in 
which  the  rate-independent  equations  are  conservation  laws. 

We  assume  that  the  static  yield  surface  associated  with  the  flow  rule  A,, 
is  given  by  Eq.  (3.70).  In  particular,  |[A||  >  0  when  <p  >  0.  so  that  we  can 
define  the  relaxation-time  parameter  r  >  0  to  be  r  —  <p/||A||.  (For  in¬ 
stance.  in  the  Steinberg-Lund  model,  r  is  proportional  to  C2.)  Noting  that 
A,-  -  =  0  and  h  —  0  when  <p  <  0,  Eqs.  (3.66)  and  (3.67)  can  be  written 


(MP)(| 

L,k 


1 

T 


elt,fky, 

K,V) 

|  +  *u(e*i'e*i'fky'K'  rj)’ 

epklJky , 

K.rj)] 

|  +h(ek„e\ (5.2) 
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where  ||A||  =  1 .  Here  [  ■  ] .  denotes  the  positive  part  of  its  argument:  [^] .  is 
(f  when  ip  >  0  and  is  zero  otherwise. 

If  the  material  is  in  the-rplastic  region  <p  >  0,  then  Eq.  (5.1)  dictates  that 
e’’  relaxes  to  the  yield  surface  <p  =  d  on  the  time  scale  of  r.  In  the  limit 
r  -*•  0,  the  How  process  is  approximately  rate-independent.  This  approxi- 
mation  is  valid  when  the  time  scale  fc  returning  the  plastic  material  to  the 
yield  surface  is  much  smaller  than  the  time  scale  for  other  How  processes. 
The  limiting  equations  can  be  formulated  as 

(L;cv),,  =  nAn(eu,ef,.fky,K,  ?? )  (5.3) 

L,k  =  nh(eu,c£,,  f\,K.  77).  (5.4) 

where  p  is  a  /.ngm/ige  multiplier  that  is  subject  to  the  loading / unloading 
conditions 


p.  >  0.  <£  <  0,  pip  —  0.  (5.5) 

written  as  a  Kuhn-  fucker  complementarity  condition,  and  the  persistence 
condition  that 


L.  f  =  0  if  p  >  0.  (5.6) 

Remark.  This  formulation  has  been  extensively  discussed  by  Simo  and 
coworkers  [30.  48.  39].  Conditions  (5.5)  can  be  derived  as  optimality 
conditions  for  a  constrained  optimization  problem:  maximize  the  rate  of 
dissipation  of  energy  during  plastic  flow,  subject  to  the  constraint  ip  <  0. 

In  the  following  we  describe  how  the  rate-independent  equations  arise 
from  the  rate-dependent  equations  as  7  -»  0.  assur  ing  that  the  How  is 
smooth.  Comparing  Eqs.  (5.1)  and  (5.2)  with  Eqs.  (5.3)  and  (5.4)  suggests 
the  identification 


1 

M  *  Hm  —[(/?]*.  (5.7) 

t  — *  0  T 

Therefore  we  calculate  lim.  r~  '[<pL. 

As  a  preliminary,  note  that  the  yield  function  ip  satisfies  the  equation 


a 


L,  ip  =  [  L,  <p  ]clasllc - [<?]  v 


(5.8) 


where 


dip  1  dip 

[  'P  LliMic  ~  7T ~(Lie)'l  ~~  77T  737  (l  ,k 


<)e. 


p6  dr, 


(5.9) 
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and 


(U f  r.  1  dtp  Hr 
iV.p  I)  Hr]  He}] 


Hip 

<>K 


1  Hi f  <<y 
t)  d'i]  dK 


h. 


(5.10) 


Equation  (5.8)  can  he  derived  by  first  calculating  T  in  Lagragian  coordi¬ 
nates  and  then  transforming  to  Eulerian  coordinates. 

Remark.  The  combinations  of  derivatives  of  c  that  appear  as  coeffi¬ 
cients  in  the  definition  of  a  are  the  derivatives  of  e  with  t .  rather  than  rj. 
held  fixed. 

At  this  point  we  must  assume  that  a  0:  this  constitutive  assumption  is 
a  stability  condition,  in  that  when  the  material  is  beyond  the  yield  surface, 
i.e..  if  >  0.  the  second  term  in  Eq.  (5.8)  drives  the  material  toward  the 
yield  surface  =  0. 

We  regard  Eq.  (5.8)  as  an  ordinary  differential  equation  for  if  along  a 
particle  path,  with  [ /.  cLi. serving  is  a  known  forcing  term.  To  simplify 
the  analysis,  we  assume  that  t  has  been  rescaled  so  that  a  is  constant. 
Furthermore,  we  suppose  that  the  flow  begins  within  the  yield  surface,  i.e.. 
if  takes  on  the  (r-independent)  initial-value  <  0  at  time  i  -  (1.  In 

order  to  solve  the  initial-value  problem,  we  consider  three  cases. 

First,  if  if !,.  .I,  =  (i  and  [/.,  if  ]obMlti, -o  >  (k  then  <f  >  0  for  small  t  >  0  so 
that  (VL=  f.  Therefore  multiplying  Eq.  (5.8)  by  expinf/r)  and  integrat¬ 
ing  both  sides  along  a  particle  path  shows  that  if  is  given  by 


*  -  f*  P-fl . «<&•  (5-U) 

-  o 

Note  that  when  t  >0  is  fixed.  (<r/r)exp[-«(/  -  v)/r]  approaches  a 
delta-function  at  .v  =  t  as  t  — *  I).  Consequently,  Eq.  (5.11)  implies  that 


1 

lim  lim  —  [g>]  _ 

r  >  T 


I 

[  ■*£  Ll.lMIC 

a 


(5.12) 


(Obviously,  taking  the  limits  t  — *  (C  and  r  -*  0  in  the  opposite  order 
yields  zero:  thus  /u.  is  discontinuous  at  t  =  0.) 

Next,  if  either  *£!,„„  <  0  or  if  =  0  and  [  L.  ^J0|-ls,lti,  <,  <  0.  then 
if  <  0  for  small  t  >  0.  Therefore 


1 

lim  lim  —  [</?]  *  =  0.  (5.13) 

, -ir  r  —II  T 

Finally,  suppose  that  if  I,  =  0  and  [  l.  ]L.|,lM1J,  .>  =  0.  Since  [<f],  >  c. 
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the  argument  leading  to  1  q.  i5.I  I)  vickis  the  inequality 

U']  «'i»\  n.  /'«■  ' [  /-  v  Ju.iMK  </s 


(  5 . 1 4  i 


Thus 


I 

ti  •  lim  imisup  [ej.  -  m;ix{().  [  /,  c  ]  .  fc  ..}  =-  t|,  (5.151 

•  U  -  .  (  i 


As  concerns  the  Lagrange  multiplier  ,u  •  0  oj  Lq.  .  we  draw  the 
following  conclusions.  Equation  (5. 1 5)  says  that  /a  (life  •  0.  Moreover, 
by  l;qs.  (5.12).  (5.15).  nad  (5  15).  p  0  only  it  Lq.  1 5 . 12)  holds,  in  which 
ease  Lq.  (5.5)  implies  (hat  /  o.  In  other  words.  Lqs.  (5.5)  and  (5.4) 
and  conditions  (5.5)  and  (5.(0  are  obtained  in  the  rate-independent  limit. 

These  conditions  lead  to  the  strain-space  formulation  of  the  loading 
criterion  for  a  rate-independent  material  [22],  We  define  p[nt]  to  be  the 
right-hand  side  of  Lq.  (5.12).  i.e.. 


M  irui 


«'t£ 

— "  (  l-  <■'),, 


1  (V 
pH  a  17 


j  ll  if 

!  a  if  af 

1 

i  H<f 

I  _ 

i  Hif  (if  i 

1  ."v7 

II  <>T]  *V,P, 

*  it 

1  <v 

tl  (1 T)  <>K  1 

(5.1b) 


and  say  that  the  How  is  undergoing  loading  when  c  =  0  and  julnil  >  0. 
neutral  loading  when  e  -  d  and  p,n<]  -  0.  unloading  when  c  =  I)  and 

/utri;il  0.  and  elastic  deformation  when  c  •  0.  Then  the  Lagrange  multi¬ 
plier  p  that  appears  in  Lqs.  (5.5)  and  (5.4)  is 

I  utril,  in  the  case  of  loading. 

^  (5.17) 

(0  otherwise. 

This  classical  formulation  is  implied  by  conditions  (5.5)  and  (5.6).  as 
follows  [29].  First  suppose  that  if  =  0  and  ulrill  >  0.  Then  in  order  for  c 
to  remain  nonpositive,  we  must  have  that  Lp  <0.  Since  L  p  = 
[L,  if  ]t.|as(ic  -  pa.  this  implies  that  p  >  /atrM,  >  0.  By  the  persistency  con¬ 
dition.  we  conclude  that  L.p  =  0.  i.e..  [Lt  <p], ,ilsIlt  -  pa  =  0.  By  the 
persistency  condition,  we  conclude  that  L,  if  =  0.  i.e.,  [L,  e]c,.lMK.  -  pa  = 
0.  Therefore  p  =  /ama,.  Next  suppose  that  <f  <  0.  Since  pup  —  0.  this 
implies  that  p  =  0.  Finally,  if  triai  <  0,  then  p  -  0  as  well:  otherwise  p 
would  be  positive,  implying  that  L,tf  =  0.  which  would  entail  that  p  - 
/i.,nal,  a  contradiction. 
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The  classical  flow  rule  is  inadequate.  however,  when  the  How  is  not 
smooth.  In  general  we  wish  to  allow  Hows  containing  discontinuities,  i.e.. 
solutions  with  boundecf-variation  that  satisfy  the  governing  equations  *n 
the  weak  sense  (the  sense  of  distributions).  In  such  a  case.  Eq.  (5.16)  might 
involve  multiplying  a  discontinuous  function  by  a  distribution  with  a 
delta-function  singularity,  a  process  that  is  ill-defined.  In  contrast,  the 
Kuhn-Tucker  conditions  (5.5)  are  sensible  for  bounded-variation  solu¬ 
tions:  the  conditions  hold  almost  everywhere.  Furthermore,  if  the  persis¬ 
tency  condition  (5.6)  is  replaced  by  the  equivalent  condition 

it 

—  (pip)  +  (pipck).k  =  0  in  the  interior  of  the  support  of  ya •  (5.18) 

fit 

then  it  also  can  be  required  of  weak  solutions. 

Thus  the  weak  formulation  of  the  governing  equations  for  rate-indepen¬ 
dent  plasticity  consists  of  the  conservation  laws  (2.22).  (2.23)  and  either 
(2.24)  and  (2.25)  or  (2.26),  together  with  the  flow  rule 


+  (p/V,p,/V‘).*  =  P/*/'A,/V  (5.HM 

d 

~{px)  +  {pKr  ) .  *  =  pfih.  (5.20) 

at 

where  /a  is  constrained  by  conditions  (5.5)  and  (5.18).  We  call  the  resulting 
system  a  constrained  system  of  conservation  laws. 


6.  Wave  Structure 

In  this  section  we  give  a  preliminary  analysis  of  scale-invariant  solutions 
for  rate-independent  plasticity.  In  the  Lagrangian  picture,  the  governing 
equations  (with  heat  conduction  neglected)  are 

f'.  = 

(6.1) 

-  (C-v'Tf 

(6.2) 

pM(\yy‘  + 

(6.3) 

= 

(6.4) 

K  =  pH. 


(6.5) 
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where  y  is  subject  to  the  constraints 

y  >-.(),  <t>  <  0.  =  0.  (6.6) 

d>  =  0  if  y  >  «.  (6.7) 

We  seek  scale-invariant  solutions  of  these  equations.  By  this  we  mean  that 
the  flow  variables  I  S.  £,^,  and  K  depend  solely  on  t  ■-  V,  A'"//, 
whereas  y  has  the  form  /!(£)//. 

If  such  a  solution  has  a  jump  discontinuity  at  speed  .S\,  then  each 
(distributional)  derivative  appearing  in  Eqs.  (6.1M6.5)  is  proportional  to 
the  delta-function  8(NyXy  -  Ss,t)  supported  on  the  plane  £  =  ,S\.  For 
instance.  E^fi  -  -Sx  A Eftti  8  and  V  (i  =  N(t  AG'  8.  The  source  terms  on 
the  right-hand  sides  of  Eqs.  (6.4)  and  (6.5).  however,  are  bounded  varia¬ 
tion  functions,  with  no  delta-function  component.  Therefore  the  jump 
conditions  (4.6)— (4. 10)  are  satisfied  (with  A Qy  =  0).  subject  to  the  con¬ 
straint  that  <i>  <  0  on  both  sides  of  the  discontinuity.  In  other  words,  jump 
discontinuities  are  unaffected  by  the  plastic  source  terms.  There  is  no  need 
to  determine  y,  in  fact,  since  the  plane  £  =  Sv  has  measure  zero,  the 
value  of  the  distribution  y  on  this  plane  is  indeterminate. 

For  a  smooth  solution,  on  the  other  hand.  Eq.  (6.2)  becomes 


PreS'  =  ~  Pr.fTF'J^S 


+  F‘ 


dSnti  <'iSa>i 

EZt.u  +  —  K 


u-p  a 

UL*yfi 


<m 


(6.8) 


Here  A!fiki>  is  the  adiabatic  elasticity  tensor  and  the  I  '0  are  the  Griineisen 
coefficients: 


A*"  ==Pref- 
_  YYa^  ■— 


d\ f 


tlF,ff  <>F^ 

d2<f 


dE„,  dS 


(6.9) 

(6.10) 


The  elasticity  tensor  is  related  to  the  adiabatic  elastic  moduli 

d2<? 

C“^5  -  pref 


^Faji  dEyb 


(6.11) 


through  the  equation  A!pks  = 
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Furthermore.  Eq.  (6.3)  can  be  replaced  by  Eq.  (3.18),  i.e.. 


96'  . 

A 


<1EP 


Jk 


H 


This  implies  that  <6  =  ‘J)|ei.1N,IC  -  na  with 


<t>U, 


iHt> 


and 


(6.12) 


(6.13) 
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/Ml> 

1  /Ml) 

1  /mi> 

1  /Ml>  fitf' 

;>  e  p 

T  TV 

- 

in-::, )  "0 

(/w  " 

T  ihS  /'  K 

(6.14) 


Therefore,  when  a  scale-invariant  solution  is  differentiable,  Eqs. 
(6. 1 )— (6.5)  become 


-e(F'„h-W)(, 

P*<4(Y'h  -^»(F,s)f  -  P^TFWS' 

dS“e  <IS°“ 

r  ( )^—k, 


+  F'„N, 


dEp- 

1  '-'vrt 


-  asc  = 


d#  .  a#  . 
A„«  +  —  H 


dEp,  '  dK 


€{£%(})( 

-  ^  =  HH. 


If  we  let 


1  d0> 


9K 


^  trial  ,= 


«^‘/y  " 


(6.15) 


(6.16) 

(6.17) 

(6.18) 

(6.19) 

(6.20) 


then  according  to  conditions  (6.6)  and  (6.7),  /I  =  /Itria|  when  <!>  =  ()  and 
Mtriai  >  0’  whereas  fx  =  0  otherwise. 

When  £  —  0,  Eq.  (6.15)  implies  that  (V')^  =  0,  while  Eqs.  (6.1 7)— (6. 19) 
are  solved  by  taking  jl  -  0,  Thus  there  is  no  velocity  variation  and  no 
plastic  flow  at  £  =  0.  Equation  (6.16)  represents  three  linear  conditions 
connecting  the  remaining  unknowns  (Fj6)€,  S(,  (Epfi)i,  and  K(.  These 
conditions  are  the  infinitesimal  version  of  the  requirement  that  the  normal 
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stress  F‘ttSi‘liNlj  be  continuous  across  a  discontinuity  propagating  at  speed 
Sy  -  0  (cf.  Eq.  (4.6)). 

On  the  other  hand,  when  £  ^  0.  Eq.  (6.16)  can  be  replaced  by 


PnttHV'h  •{As)\{V>)e  +  fiFlaNp 


dS 
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i>Ey« 


US  "  * 

7F 


+  /v.,r^ 


7  \ 


^  Prel 


f-a/J  . 


Jk 


H 


(6.21' 


where  A,;v  =  /V^A'^A^  is  the  acoustic  tensor.  Since  /i,nii  is  proportional  to 
( )f ,  Eq.  (6.21)  represents  an  eigenvalue  problem:  if  <I>  =  0  and  /2tn:il  >  0. 
then  prcl £ 2  is  an  eigenvalue  of  the  plastic  acoustic  tensor  (\vx  Yl  defined  by 


J  a 


I  dSali  tic? 

+  prcfr^ 


;)EU  / 


A. 


+ 


dSafi  d<?y 

IF  +PrefP^, 


H 


BE 


■FI(NV:  (6.22) 


ev 


and  otherwise  pre[£2  is  an  eigenvalue  of  the  (elastic)  acoustic  tensor 
(Ay)' .  The  eigenvector  in  each  case  is  ( V')i ,  and  once  it  is  known,  the 
remaining  unknown  quantities  (Pa)f,  Sf,  (E^p)f,  and  K ^  are  determined 
by  Eqs.  (6.15)  and  (6. 1 7)— (6. 19). 

Assuming  the  eigenvalues  of  (AiV)';  and  (A^)',  to  be  positive,  the  wave 
modes  occur  in  three  pairs.  For  instance,  in  the  elastic  flow  of  an  isotropic 
material,  the  fastest  wave  corresponds  to  longitudinal  stress  (i.e„  pressure) 
modes,  while  the  two  slower  waves  correspond  to  radial  shear  stress  (i.e.. 
necking)  and  angular  shear  stress  (i.e.,  torsion)  modes.  In  plastic  flow 
(p  >  0),  waves  propagate  at  speeds  corresponding  to  the  eigenvalues  of 
the  matrix  (A^)',,  which  is  a  rank-one  perturbation  of  (A  v)';.  Typically  the 
constitutive  model  leads  to  plastic  wave  speeds  that  are  smaller  than  the 
corresponding  elastic  wave  speeds. 


7.  Conclusions 

In  this  paper  we  have  formulated  the  equations  governing  rate-depen¬ 
dent  and  rate-independent  plastic  flow  of  metals  in  a  conservative  form, 
applicable  to  both  smooth  and  discontinuous  solutions.  A  preliminary' 
analysis  of  the  wave  structures  arising  in  the  rate-independent  case  was 
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also  given.  This  work  provides  the  foundation  for  further  studies  of 
nonlinear  waves  in  plastic  flow,  and  for  the  use  of  this  information  in 
numerical  computations. 
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